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Abstract 

<n: 

t5 We describe how thermalization occurs in heavy ion collisions in the frame- 

work of perturbative QCD. When the saturation scale Q s is large compared 
to Aqcd, thermalization takes place during a time of order a" 13 ' 5 ^^ 1 an d 
the maximal temperature achieved is a 2 ' 5 Q s - 

m ■ 

>: 

I. INTRODUCTION 

It is possible that at RHIC, for the first time, heavy ion collisions occur at energies 

high enough to be described by perturbative QCD. At the Large Hadron Collider (LHC) 

perturbative QCD is expected to work even better. At these energies we will assume that 

^ ■ immediately after the collision the initial distribution of gluons is given by the saturation 

scenario |lp||. Thus the relevant hard scale is the saturation scale Q s , estimated to be 1 

gn! GeV at RHIC and 2-3 GeV at LHC. 

The single most important question in the physics of heavy ion collisions is thermaliza- 
tion. The conventional argument in favor of thermalization is that at larger collision energy, 
more gluons are freed in the first moment after the collision, and these gluons collide more 
frequently with each other. However, the distribution of these gluons is initially very far 
from thermal equilibrium. In addition, the strong coupling constant decreases at high en- 
ergies, making thermalization harder to achieve. Whether the system has enough time to 
equilibrate before falling apart is thus a delicate question requiring detailed consideration of 
different physical processes. 

In this paper we show that, in the limit Q s ^> Aqcd corresponding to very large nuclei 
and/or very high collision energy, thermalization occurs relatively fast while the system is 
still undergoing one-dimensional expansion. The unexpected feature of our analysis is the 
way thermalization occurs. During the first period of time the most important process is 
the emission of soft gluons which overwhelm, in terms of number, the primary hard gluons 
at time r ~ a"^' 2 Q~ l . These soft gluons then quickly equilibrate and form a thermal bath, 
which initially carries only a small fraction of the total energy. The thermal bath then 
draws energy from the hard gluons. Full thermalization is achieved when the primary hard 
gluons have lost all their energy. Parametrically, this happens at r ~ a^ 13 ' 5 Q^ x , at which 
time the temperature of the system achieves the maximal value of a 2 ^Q s . Surprisingly, the 
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time dependence of the temperature of the soft sector in the (admittedly narrow) region 
0,-5/2 ^ q^ t ^ a - 13 / 5 can be found almost analytically, and the result depends only on 
the total number of primary hard gluons. It is important to emphasize here that we only 
provide parametric qualitative and quantitative estimates based on the assumption that 
always a < 1. A more realistic study is required in order to figure out the numerical 
coefficients of the given estimates. 

This picture of "bottom up" thermalization is different from that considered in a previous 
work by one of us (see also 0) which did not take into account particle production, and 
also of Refs. |@-H|. In [1 Ij the importance of inelastic processes, even for kinetic equilibra- 



tion, has been observed when using the relaxation time approximation for the collision term 
and with initial conditions different from the ones considered here. Although the analysis 
requires a small coupling, one can hope that many qualitative features of this picture are 
already present in heavy ion collisions at LHC or even RHIC energies. At the very least, 
the finding gives us confidence that thermalization always occurs in heavy ion collisions at 
sufficiently high energies. Our main emphasis here is the thermalization that begins with the 
softer momentum modes. Still more detailed investigations, including more specific realistic 
initial conditions, are required for treating the high momentum tails of the distributions 
which are expected to require a longer thermalization time, but which are important for 
plasma signatures sensitive to momentum scales larger than Q s . 

II. QUALITATIVE DESCRIPTION OF EARLY TIMES 

We will be interested only in the central rapidity region of central collisions. In this region 
one can assume boost invariance: all physical quantities depend only on the proper time 
t = a/£ 2 — z 2 , but not on the rapidity rj = |ln |^, where z is the direction of collision. For 
large nuclei, the dependence on the transverse (x and y) coordinates can also be neglected, 
which is equivalent to assuming the medium to be infinite in the transverse directions. With 
these simplifications, all physical quantities depend only on the single coordinate r. The 
evolution proceeds through several regimes, with Q s r ~ a -3 / 2 , a _5//2 , and a~ 13 / 5 marking 
the borders between neighboring time periods. 

We will describe qualitatively the evolution of the system up to Q s r ~ a" 5 / 2 . In all 
subsequent estimates, we will assume In - ~ 1 and keep track only of powers of a itself. 
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A. Very early time, 1 <C Q s t <C a 



-3/2 



At the earliest time, r ~ Q^ 1 , gluons are freed from the nuclei. These gluons have typical 
momentum of order Q s and occupation number of order 1/a. This is the basic assumption 
of the rest of this paper. Due to the large occupation number these gluons interact so 
strongly that it is more appropriate to describe them as a nonlinear gluon field rather than 
a collection of particles. Only when Q s t becomes larger than 1, the classical field becomes 
almost linear, and one can start to describe the gluons as particles on mass shell with a 
well-defined distribution, 



J*_ = l f (kL) (1) 

dydkl a\Qj' K ' 

The precise form of / depends on the details of the nonperturbative physics at the scale 
t ~ Qj 1 and is not a subject of the present paper. A promising method to find / is by 
simulation of the classical gluon field [|1| . We will call the particles produced during these 
first periods the hard gluons, since later gluons with smaller momenta will be produced. 

The density of hard gluons [|l]-|6|] decreases with time due to the one-dimensional expan- 
sion, 

**~-7^v (2) 

Oi{Q s T) 

If there was no interaction, the occupation number would remain of order 1/a due to the 
Liouville theorem. There is no contradiction with Eq. (Q), since the typical longitudinal 
momentum of hard gluons also becomes smaller (p z ~ 1/r) because gluons with larger 
longitudinal momentum escape from the spatial region under consideration during a time of 
order r. 

In reality, gluons interact by elastic scattering. Most of the scatterings are small angle, 
with exchange momentum q <C Q s . The effect is the broadening of the distribution along the 
p z direction, thus lowering the typical occupation number. The lowest possible momentum 
exchange is the Debye mass, which is determined by fll6|,|8|,[ri 



2 / J3 



mn ~ a 



f h {p) aN h Q 



d 3 p^~^~-^. (3) 

v Qs Q s t 



If one assumes that mo "C p z (this condition will be verified a posteriori) , most collisions 
do not take particles away from the momentum region where the occupation number is 
large (k± ~ Q s , k z ~ p z ). The frequency of collisions that a typical particle encounters is 
enhanced by the Bose factor, 

c^V co i aN h 

— aN h (l + f h ) 2 > 4 

dr mf)p z r 

where a ~ a 2 m^ 2 is the cross section, and fh = Nh/{Q 2 s p z ) is the typical occupation number, 
which is assumed to be large. These random collisions increase the longitudinal momentum 
of gluons, which is typically 

- 2 - ■ AT col m 5 ~ ^ , (5) 
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which implies 



P-W'-^. (6) 



This relation may be derived from 0, taking the Bose enhancement into account, which 
amounts to replace the average p\ ~ aQ 2 s by p 2 ~ a s Q 2 {l + fh) ~ aN h /p z . ^From Eqs. 



(H) and (|>D we see that mo <C p z , as assumed. The typical occupation number fh ~ 
®~ 1 {Q s tY 2 ^ is large until Q s t ~ a -3 / 2 . 

Beside elastic scatterings there are also inelastic scatterings in which gluons are produced. 
As we will see, the most important produced gluons are those with smallest energies. In 
principle gluons with energy as low as mo can be produced. However, once produced, the 
momentum of these gluons is pushed up by multiple elastic scatterings with hard gluons. 
As a consequence of (|B|) the smallest momentum of soft gluons is of order p z : k s ~ p z ~ 
Qs/(Qs T Y^ 3 - The number of soft gluons with momentum k s at time r that are produced at 
this moment is estimated from the Bethe-Heitler formula to be (|TS| 



N , „ T ^ . rftiM^l + A) 2 ~ r^m + f„f ~ ^£^ . (7) 

The time interval around the moment r is taken to be of order r. Once these soft gluons 
are produced they will remain, and their density is decreasing as 1/r. One notices that 
N s /k s ~ Nh/Q s , so the Debye mass receives equal contributions from hard and soft gluons, 



., aN h x aN s 



m ^ (8) 

but the estimate of Eq. (||) is still valid parametrically. 

B. Setting up the stage for thermalization: a 3 / 2 <C Q s t <C a 5 / 2 

Beginning from Q s t ~ a -3 / 2 , the occupation number of hard gluons drops below 1, and 
the estimates of the previous section need to be revised. We will see that when a~ 3 / 2 <C 
Q s t "C a -5 / 2 soft gluons contribute negligibly to the total number of gluons but give most 
of the Debye screening. In other words, 



N s < N h , 
ml ~ -j— • (9) 

First let us estimate k s , which is the typical momentum of soft gluons. When fh <C 1, we 
have 

k\ ~ N col ml ~ raN h m 2 D ~ aQ 2 s , (10) 

which is now a constant scale. The number of soft gluons that have been produced at time 
t is 

N g ~T^N*~^. (11) 



^From Eqs. (m and (11), one finds 
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m D 



a 3 
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For the soft gluons to give the dominant contribution to Debye screening one needs N s /k s ^> 
Nh/Qs, which requires Q s t ^$> a -3 / 2 . The number of soft gluons become comparable to that 
of hard ones at Q s t ~ a~ 5 ^ 2 . Therefore all estimates in this section are valid in the interval 
o^/ 2 < Q s r < a- 5 / 2 . 

In-medium emission of gluons is suppressed by the Landau- Pomeranchuk-Migdal (LPM) 
effect. The latter is operative at scales larger than A; L pm determined by k LPM = m^iV/jer) -1 



or 



fcLPM ~ ^ ~ <* 1/2 Qs > ( 12 ) 

a 2 N h 

which is parametrically of the same order as the scale k s . The LPM effect suppresses the 
gluon production compared to the Bethe-Heitler rate at scales larger than k s , but does not 
affect the rate of producing particles in the interval between run and k s . Therefore all 
previous estimates based on the Bethe-Heitler formula remain intact. 

III. THERMALIZATION OF THE SOFT SECTOR: Q s t > a" 5 / 2 
A. Qualitative description 

After Q s r ~ oT^I 2 most gluons are soft, N s 3> Nh- We will see that the soft gluons collide 
very frequently with each other, and achieve thermal equilibration amongst themselves. The 
soft sector is characterized by the temperature T, which is a function of time. The system 
as a whole is still not in thermal equilibrium, since most of the energy is carried by a small 
number of hard gluons. These few gluons collide with the soft gluons of thermal bath and 
constantly loose energy to the latter. 

A hard gluon with energy of order Q s looses energy to the bath by the following mecha- 
nism. First it emits a particle with a softer momentum fcbr, which, during a time comparable 
to r, splits into two gluons with comparable momenta (hard branching). The products of 
this branching quickly cascade further, giving all their energy to the thermal bath. 

As we will verify, k^ r lies in the region where the emission rate is LPM-suppressed. One 
can estimate the time of emission, tb r , for a gluon having momentum k^ as follows: 1/tbr ~ 
a/tf, with the formation time U given by if ~ k^/k 2 . The gluon picks up a momentum k t , 
transverse to its direction of motion, given by k\ ~ rri^U/X. Using A -1 ~ N s a ~ N s a 2 /m^ 
one obtains (see also Eq. (|T9|) below), 

1 a 2 N 1 / 2 



Equating the branching time tbr with r, and using N s ~ T 3 , where T is the temperature of 
the soft thermal bath, we find that fcbr ~ a 4 T 3 r 2 . The time dependence of the temperature 
T is still to be found. 

The number of &br-gluons produced per unit time per unit volume is 

dN(k hT ) N h a 2 N}/ 2 N h Q 2 

dr t hl k ]l 2 ar 2 ' { } 

Subsequently, the rate of energy flow from the hard gluons to the soft thermal bath is 

fc br — -; ~ a Q s l . (15) 

This energy flow increases the energy in the thermal bath, and thus must be proportional 
to d(T 4 )/dr. Therefore one finds 

T~a 3 Q 2 s r. (16) 

The temperature of the soft thermal bath increases linearly with time, even when the system 
is expanding, due to the hard gluons which serve as an energy source. When the bath 
starts forming (Q s t ~ a -5 / 2 ), its temperature is a l / 2 Q s . The relaxation time of the soft 
sector is of order r re i ~ (a 2 T) _1 ~ (a 5 Q 2 r) _1 . When Q s t > a~ 5 / 2 , r rc i -^ r, which 
justifies the assumption of thermal equilibration of the soft sector, i.e. in the relaxation time 
approximation f s — > / eq (l — exp(— r/r re i)). 

The linear growth of T terminates when the hard gluons loose all of their energy. This 
happens when k^ T ~ Q s , or r ~ a -13 / 5 ^^ 1 , when the temperature achieves a maximal value 
of order a 2 ' 5 Q s , which is larger than the initial temperature only by a factor of a' 1 ' 10 . 
Subsequently the temperature decreases as r -1 / 3 [^ . 

The border between the Bethe-Heitler and LPM regimes is at /clpm = ^d/^ 2 -^^ ~ 
T ~ a 3 Q 2 r. On the other hand k^ T ~ a A T 3 r 2 ~ a 13 Q®r 5 , which is much larger than /clpm 
when Q s r ^> or h l 2 . Thus the use of the LPM formula in Eq. (O) is justified. 



B. Quantitative description 

The picture given above can made quantitative in the form of a Boltzmann equation 
describing the kinetics of the branching process. In this kinetic approach the evolution of 
hard modes follows the equation: 



d p z d 



,<9t t dp 
where C e \ is the elastic collision integral (see Ref. J/J) and 



)f(p) = C cl + C pvod , (17) 



1 d 2 I 



r ^= i< L <— ^(4ik(-)( i +/(^)( i +/(^ L ^))-/®/(^ L ^)( i +/(-)) 

; lx 0/z ix x xx 



yprod " J X dxdtXx^ 2 



1 



/(p)(l+/(px))(l+/(pU-s))) " f(px)f(p(l-x))(l+f(p)) 



(18) 



is the term describing the 2 — » 3 and 3 — > 2 processes. The four terms in the collision integral 
correspond to the diagrams (a), (c), (b) and (d) in Fig. [l]. The black blob in this Figure 
represents multiple scatterings off individual gluons in the medium. In Eq. ([18]) d 2 I/dxdt is 
the rate of a hard gluon with momentum p ~ Q s to split almost collinearly into two gluons 
with momenta px and p(l — x) while moving in a medium where most particles carry much 
smaller momenta. This rate is LPM-suppressed and can be computed using the method of 
Ref. [|T9"|-|2"T|] , for an infinite size medium, which yields 



where 



h(x) = h 



and 



cPl 

dxdt 



a 2 iV 1 /2 



P 



1/2 



h(x) 



[l-x + x 2 fl 2 
[x — x 2 ) 3 / 2 



N 2 



TT 1 / 2 (AT2 -1)1/2 



hi 



<*?> 



m 



D 



1/2 



(19) 



(20) 



N = 2(N 2 - 1) 



(2tt) 



:/(p)(l + /(p)). 



(21) 



The A^ in (21) represents the density of possible scatterers in the system. N is not quite the 
same as the density of gluons because of the extra factor of (1 + /), a factor necessary to 
correctly give the interaction rate. From the discussion given earlier we expect N s ^> Nh so 
long as Q s t > a~ 5 ^ 2 and so we expect iV to be dominated by the soft (thermalized) particles 
in our present discussion. 

h(x) has a symmetry property that h(x) = h{l — x), using which one can check that 
the collision term (p~8|) conserves energy, JdppC pro d(p) = 0. One might wonder why C pro d is 
of order a 2 instead of a 3 . An intuitive derivation for this dependence is already given just 
above Eq. (pi). 
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FIG. 1. The diagrammatic representation of C pro d 



A remark on the Boltzmann equation is in order. In Ref. a Boltzmann equation 
without particle production was considered. The system does seem to approach kinetic 
equilibration [0|§, but during a relatively long time (~ exp(a _1 / 2 )Q7 1 ). This is because most 
of the elastic scatterings are at small angle, and it takes a lot of small angle scatterings to 
change the particle distribution considerably. In more technical terms, there is a cancellation 
between gain and loss terms in the elastic collision integral. The 2 — > 3 process, on the other 
hand, is asymmetric and does not have this type of cancellation, and so is important although 
it is smaller than the elastic rate. Parametrically higher-order processes like 2 — > 4 can be 
neglected since they are suppressed by an additional factor of a compared to 2 — ■> 3 while not 
adding any qualitatively new feature to the evolution. Whether in a realistic circumstance, 
say at RHIC, higher order processes are negligible or not is beyond the scope of the present 
discussion |IU |. 



The elastic collision integral C e \ is responsible for the thermalization of the soft sector 
p ~ T, and for p z broadening of the hard particles. The inelastic term C pro d is responsible 
for the energy flow from the hard to the soft sector. In principle, the products of branching 
of hard gluons have a small transverse momentum, but this is negligible compared to the 
subsequent broadening by elastic scattering. Therefore, C pro d is written as if the branching 
was exactly collinear. 

For hard gluons which have small occupation number (/ <C 1) Eq. ([D]) simplifies con- 
siderably. If we introduce 



e(p ± ) = 27cpljdp z f(p) (22) 

which is normalized so that the energy carried by hard gluons is Jdp±e(p±), then the kinetic 
equation for e(p±) becomes 



1 d a 2 N l > 2 , 

-^( re (^)) = — 172 - ldxh(x) 



^JPA_l 



* 1/2 * — -^(p. 



x 



(23) 



In the following we will write p instead of p± for simplicity. Eq. (^) does not describe the 
soft gluons, for which C e \ is essential. Fortunately, these gluons are fully equilibrated and 
one can characterize the soft sector by the temperature T. In particular, 

N = g N T 3 , g N = 2{N 2 c -l)\. (24) 

The time evolution of T depends on the amount of energy that flows from the hard to the 
soft sector. To find this quantity, let us introduce an intermediate scale p , T <C p <C Q s , 
and integrate Eq. (p3[) from p to 00. One finds 

„ OO 1 Po/x 

-^-(rldpeipyj = -a 2 N 1/2 fdxh(x)x J dp^. (25) 

po Po 

The left hand side of Eq. (|25| ) has the meaning of the rate of energy flow from above po 
to below p . This quantity must be a constant independent of p Q when p Q <^C Q s . This is 
possible only when e{p) has the following behavior at small p: 

8 



e(p) = -^2, T<€.p<€.Q s . (26) 



The right hand side of Eq. (|25| ) is then — 6/i a 2 A^ 1 / 2 e 1 , where 



/■ fl _ x + x 2 \5/2 1 

6= /(fa iL g + *J g in -^4.96. (27) 

v ' 



The equation that governs the evolution of the temperature is then 

1 d 



(r 4 / 3 e s (T)) = bh^N 1 ' 2 ^ (28) 



because of the one-dimensional expansion [0,0 where e s (T) is the energy density of the 
soft gluons at temperature T, 

e s (T)=g E T\ g E = 2(N 2 - 1)^ . (29) 

Eqs. (0), (p4|), (|28|), and (0) are the equations that govern the evolution after Q s r ~ a -5 / 2 . 
From Eq. (^) one needs to extract ei, which enters Eq. (^g) that describes the temperature 
T, which feeds back to Eq. (0) through iV. 

Let us show that this procedure can be performed analytically in the regime a -5 / 2 <C 
Q s t <C a -13 / 5 . Eq. (0) can be solved by iteration. We will be looking for the solution in 
the form 

e(r,p) = eo(r,p) + £i(r,p) + e 2 (r,p) H (30) 

where €q(t, p) is the energy distribution which starts the iteration. It is related to the number 
of hard gluons by 

N h (r)= / dkk- l e Q {r,k) (31) 

Jp 

so long as p ^> T and p <C Q s . Due to the one-dimensional expansion it decreases as 1/r. 
Solving for ei(r,p), we find, for p <C Qs, 

ei(T,p) = ^# (32) 



1/2 



P 

where ei(r) satisfies the equation 

^-(re 1 (r))=a 2 h N 1 / 2 N h r. (33) 

We see that ex(r,p) is singular at small p, and hence at very small p, e\ ^> e . One could 
expect that at these small p's the series (0) is ill behaved. However, if one tries to do the 
second iteration and inserts (|32|) into the collision integral of Eq. (|2~3|) , the result is zero. 
This is because e{p) = const jp x l 2 is a formal static solution to the Boltzmann equation 

9 



(f23|). This solution describes a stationary state with a constant energy flow from high to low 
momenta. As a consequence, at small p the iteration procedure stops at the first iteration, 
and the solution to the Boltzmann equation is simply eo + e\. Notice that Eq. fl3"2| ) has the 
form expected in Eq. (p6|) . 

To solve Eqs. ( |33"D and ( |28| ) one needs to know only the total number of the hard gluons, 
Nh, but not the whole distribution function. This number can be parametrized by Q s and 
a dimensionless constant c 0: 

N 2 - 1 O 3 

N " = T^^ry (34) 

Eqs. (|33| ) and fl28"|) then can be solved. For the temperature, we have T = CT<y 3 Q 2 r, where 

3 N?-lg Nll2 36 „ , (A; 2 ) 

c T = -—2- — ^-6/i 2 c = — jv 3 cln^/- . 35 

160tt 2 N c g E ° 8tt 5 c ml K ' 

Eq. ( J35|) is valid only with logarithmic accuracy. This is due to the fact that the rate of 
branching by a hard particle in a thermal medium d 2 I /dxdt is known only with this accuracy. 
Since k\ ~ a 2 N s kb T , k 2 /ml ~ a 5 (Q s r) 2 , and the argument of the logarithm in Eq. ( [35D is 
parametrically large for Q s r ^$> a~ 5//2 . An improved calculation of d 2 I /dxdt would give a 
more accurate estimate for T by precisely determining the constant in the logarithm in Eq. 
(P5j). Thus we arrive at the same results, with a more accurate determination of T, as in 
part A of this section. 
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